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We show that finite size scaling techniques can be employed to study the glass transition. Our 
results follow from the postulate of a diverging correlation length at the glass transition whose 
physical manifestation is the presence of dynamical heterogeneities. We introduce a parameter 
B(T, L) whose temperature, T, and system size, L, dependences permit a precise location of the glass 
transition. We discuss the finite size scaling behaviour of a diverging susceptibility \{L,T). These 
new techniques are successfully used to study two lattice models. The analysis straightforwardly 
applies to any glass-forming system. 
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Progress in the study of the 'glass transition' has been 
slow for two main reasons. First, almost by definition, 
the dynamics of a liquid supercooled through its melting 
transition towards its glass phase becomes so slow that 
it is impossible to actually cross a possible phase tran- 
sition while staying at equilibrium. One relies therefore 
on thermodynamic or kinetic extrapolations which stop 
at the glass temperature T g where the system no longer 
equilibrates. As a consequence, the mere existence of 
a genuine transition is still questioned, even though ex- 
periments now probe more than a decade of decades of 
relaxation times. 

Secondly, the nature of the putative transition remains 
largely unknown as compared to more standard phase 
transitions. Basically, the situation is such that the for- 
mulation of a Ginzburg-Landau type of model is still im- 
possible, because one does not know what would be the 
correct order parameter, its symmetries, dynamics, cor- 
relations, etc. Similarly, at the beginning of the 90's, 
two papers apparently settled the related question of the 
possible existence of a diverging correlation length at the 
glass transition, reporting the absence of such a length 
scale 0,0. As a regrettable corollary, modern devel- 
opments of statistical mechanics such as renormalization 
group, finite size scaling, or universality classes became 
apparently useless to study the glass transition. 

In this paper, we show that finite size scaling tech- 
niques first derived from renormalization group concepts 
in the context of continuous phase transitions are use- 
ful to study the glass transition, at variance with the 
common belief. This should allow much more powerful 
extrapolations to locate a possible transition and a bet- 
ter characterization of its nature, thus making progress 
on the two important points mentioned above. 

The paper is organized as follows. We first describe 
the observations underlying our work, before describing 
the principles of the analysis. As a first application of 
the method, we then consider two well-known models 
proposed in the context of the glass transition, namely 
the one-dimensional Fredrickson- Andersen model Q and 
the three-dimensional Kob- Andersen model |4|. 



Our work is based upon the recent burst of activity re- 
lated to the definition, observation and characterization 
of dynamical heterogeneities in the dynamics of super- 
cooled liquids 0. The underlying physics is that parti- 
cles move in a cooperative manner, the more so the nearer 
the glass transition. This picture immediately suggests 
the existence of a correlation length £ (T, t) associated to 
these dynamical domains, reflecting the length scale on 
which the particles dynamics have been correlated in a 
time interval t at temperature T. We define a 'coherence 
length' t(T) = £(T,t le \), where t Ie \ = t le i(T) is the mean 
relaxation time of the liquid. Whether or not £(T) is con- 
nected to some yet undiscovered structural length scale is 
an important open problem. Crucially, however, progress 
can be made if one postulates that £(T) actually diverges 
at some — possibly zero — temperature T c which will thus 
non-ambiguously define the glass transition temperature 
one seeks to determine. It should coincide with the tem- 
perature at which t rc \(T) also diverges. 

A measure of £(T, t) in a liquid can be obtained via the 
spatial decay of the following correlation function: 

C(r,t) = (F(r + r ,t)F(r ,t)) - (F(r + r ,t))(F(r ,t)), 

(1) 

where F(r,t) — p(r, t + to)p(r, to); p(r,t) is the density at 
time t and position r, and to are arbitrary position and 
time, respectively. Alternatively, F(r,t) can be replaced 
by any two-time function, F(r,t) = A(r,t + to)B(r,to), 
where A and B are physical observables. In essence, 
Eq. is a two-time, two-point correlation function 
which measures correlations in trajectory space 0]. 
Note that the commonly discussed 'non-Gaussian param- 
eter' is related to Eq. when A and B refer to the par- 
ticle positions, but for r — 0. As such, it contains no 
information about £(T). In fact, Eq. Q captures most 
of the measurements of dynamical heterogeneities which 
amounts to first distinguish between 'fast' and 'slow' par- 
ticles and then to check how these are spatially corre- 
lated 0. The correlators F(r,t) precisely tell us how 
fast a particle is, while Eq. Q measure the correspond- 
ing spatial correlations. The use of the dynamic func- 
tion F(r, t) makes an arbitrary criterion to distinguish 
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between particles unnecessary. Variations of the correla- 
tor ||TJ have been studied. In Refs. [HQ' the density fluc- 
tuations were replaced by an arbitrarily defined overlap 
between configurations, and by particle displacements in 

Refs. HE EI- These works show that f( T >*) is indeed 
maximum for t w i re i(T), which justifies our definition of 
£(T). Also, £(T) decreases with T, although simulations 
are yet too limited to confirm a possible divergence at 
low temperature. 

The observation of a growing coherence length scale 
has deep consequences which we start to explore in this 
paper. This indeed suggests that local dynamical func- 
tions F(r, t) become long-range correlated when the glass 
phase is approached. As such they play a role similar 
to that of the order parameter in a standard continuous 
phase transition. This simple remark is central in our 
approach. Having identified a quantity analogous to an 
order parameter, it becomes a simple task to extend the 
tools developed for conventional transitions to study the 
problem of the glass transition, starting here with finite 
size scaling techniques. 

Define first the parameter B(T,L) as 
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with if = ip(t rcl (T)) = L~ d J d d rF(r,t Icl (T)), where d is 
the space dimensionality. The crucial information is that 
the average (■ • -)l is performed in a system of finite linear 
size L. By construction, B(T, L) = when the fluctua- 
tions of ip are Gaussian, since it is built from the fourth 
cumulant of the probability distribution function of ip, 
Pl{<p)- This distribution is Gaussian when L 3> £(T) 
as a consequence of the central limit theorem. How- 
ever, Pl{}P) becomes non-Gaussian when L < £(T), and 
B(T, L) is therefore non-zero in that case. Moreover, it is 
natural to expect that B(T, L) becomes a scaling function 
of the variable £{T)/L in the regime £(T), L ^> a, where 
a is the microscopic length scale of the problem, e.g. the 
particle size in a simple liquid. Also, Pl(v>) should be 
independent of L and £(T) in the regime a <C L <C £(T). 
The latter assumption implies in particular that B(T C , L) 
is a constant 'universal' number since L <C £(T C ) — oo is 
automatically satisfied at the transition. Therefore, the 
curves B(T,L) versus T for various L cross at T c . The 
definition @ is naturally inspired by the classical pa- 
per |l2j . where a numerical finite size scaling analysis 
of the Ising model was performed. In that case, <p was 
the order parameter of the paramagnetic-ferromagnetic 
transition, i.e. the magnetization density. 
We also define the susceptibility 



X(T,L) 



1/ 
T 



(<p)l 



(3) 



again explicitly retaining its L dependence. Similarly, 
one expects this function to exhibit a scaling behaviour 
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-4 -2 2 4 6 8 10 
log(f(T)/Z) = l/T-logZ 

FIG. 1: Inset: The parameter B(T,L) as a function of the 
inverse temperature 1/T for various system sizes L = 2, 4, 
■ • ■, 1024 (top to bottom) cross at the glass transition T c — 0. 
Main: The parameter B(T,L) as a function of the scaling 
variable log(l{T)/L) with t{T) = exp(I/T). Both figures are 
for the one-dimensional spin facilitated Ising model. 



for large L and t(T). When the correlator C(r, t) is itself 
a simple function of the ratio r/£(T), one gets 



X (L,T)^L d x 



(4) 



where the scaling function x(x) behaves as x(x C 1) « 
x d and x(x ^> 1) rj const. The assumption made on the 
correlator is sup ported by all known simulation results 
so far 0, 0, 0| . It is satisfied also in the two models 
studied below. Using with caution the language of crit- 
ical phenomena, this would mean that the anomalous 
exponent 77 is 77 = 2 — d, implying a simple Fisher's law 
7 = vd between the exponents of the susceptibility and 
the correlation length. This point certainly needs further 
investigation in realistic supercooled liquids. 

The suggested analysis requires extrapolations towards 
the glass phase, as all equilibrium measurements in super- 
cooled liquids do. However, the new quantities B(T,L) 
and x(^) L) defined above are much richer quantities 
than, say, an average relaxation time since they contain 
information on the full distribution P^ip). This will be 
clearly demonstrated in the rest of the paper. 

We now use the analysis suggested above to study 
two lattice models. We start with the one-dimensional 
version of the spin facilitated Ising model introduced in 
Ref. 0. Spin facilitated models have been introduced as 
caricatures of real supercooled liquids, in the sense that 
they are non-disordered spin models displaying glassy dy- 
namics with only very simple static correlations. We have 
chosen this version of the model since it is reasonably well 
characterized Moreover, as a one-dimensional 

lattice model, it is relatively easy to simulate on very 
large time and length scales, nicely confirming the ex- 
pected scaling behaviour. These results thus constitutes 
a very interesting first application of the methods de- 
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definition (j2J is the persistence of the spin s, at site 
i. The mean persistence is given by (F(i,t)} l=oo = 



exp 
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FIG. 2: Inset: The susceptibility x(T L) as a function of 
the inverse temperature 1/T for the same sizes as Fig. 
increasing from bottom to top. Main: The rescaled sus- 
ceptibility x(T, L)L~ X as a function of the scaling variable 
\og(l{T)/L) with £{T) = exp(l/T). Both figures are for the 
one-dimensional spin facilitated Ising model. 



scribed above. The model is defined by the Hamiltonian 
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where Sj = 0, 1 are two-state observables located at the 
sites of a chain of size N with periodic boundary condi- 
tions. Glassiness originates from the chosen dynamics, 
since the static properties of the non-interacting Hamil- 
tonian JSJ are completely trivial. In a standard Metropo- 
lis algorithm, the transition rates are w(si —* 1 — s,) = 
min(l, e _/3AB ), where AE is the energy change in the 
transition Sj — > 1 — s,. Here, the rates are given in- 
stead by w'(AE) — (sj_i + s i+ i)w(AE) meaning that a 
spin can flip only if it has at least one nearest neighbor 
whose value is 1. This kinetic constraint implies glassi- 
ness. From (J5J, one has (sj) = (1 + e x ' T )~ x , so that 
the density of spin 1 becomes small at low T therefore 
allowing less and less transitions. This leads to an Arrhe- 
nius relaxation time, t TB \ = exp(3/T), which makes the 
model 'strong' according to the classification commonly 
used in the glass literature. Also, the dynamics becomes 
more and more cooperative when approaching the glass 
transition at T c = 0, with a diverging coherence length 
£(T) = exp(l/T) 

We simulate this model using a continuous time algo- 
rithm in the temperature interval T £ [0.125,1.0], ap- 
proximately corresponding to 13 decades in relaxation 
times. For each temperature, a large system of N = 2-10 4 
spins is studied, and averages are performed in subsys- 
tems of finite size L £ [2, 1024] for a number of indepen- 
dent configurations ranging from 200 to 20000. Fortu- 
nately, more statistics is required at higher temperatures 
where deviations from Gaussian behaviour are smaller. 
A natural choice for the dynamic function F(i,t) in the 
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, which defines the relaxation time. 

Our results for this model are presented in Figs.n an d 
13 The inset of Fig. ^ shows the variation of B(T,L) 
as a function of 1/T for various sizes L. As expected, 
B(T, L) crosses over from B(T, L) = at high T to a 
non-zero value at low T, the locus of the crossover being 
L-dependent. The various curves clearly cross at T c = 
only. One observes in the main figure an excellent col- 
lapse of the data, showing that B(T, L) satisfies the ex- 
pected scaling form B(T, L) w B(l{T)/L), Since P L {tp) 
is bimodal when L <C £(T), it is a simple task to com- 
pute B(T C = 0,L) = B{x -> oo) = 0.566.... The inset 
of Fig. |3 presents the susceptibility x(T,L) as a function 
of 1/T for various sizes L. Again, one clearly sees the 
expected behaviours: x(T, oo) diverges when T — > as 
e 1 ^ = £(T), while for finite L the divergence is smeared 
when £(T) > L. The main figure shows the perfect col- 
lapse of the data suggested by Eq. @ when x(T, T)L _1 
is plotted as a function ai£(T)/L (remember that d = 1). 

As a second application of the techniques described 
above, we study the kinetically constrained lattice gas in- 
troduced in Ref. Q ■ The model consists of hard spheres 
on a cubic lattice of linear size N. In this case, tem- 
perature plays no role, and the relevant control param- 
eter is the density of particles, p. This lattice gas is 
also complemented by kinetic rules which make the dy- 
namics glassy. In a non-constrained lattice gas, particles 
move to an empty nearest neighbour site with unit rate. 
In the present model, particles can move to an empty 
nearest neighbour, provided the particle has fewer than 
4 neighbours before and after the move. This kinetic 
rule aims at reproducing dynamics of a supercooled liq- 
uid, where particles can hardly escape the cage formed 
by their neighbours. As a simple lattice model for the 
glass transition, it has been much studied since its intro- 
duction. Simulations have reported a dynamical arrest 
at a density p c ks 0.881 where the relaxation time appar- 
ently diverges as a power law Q. Very recently, however, 
a lower bound for the diffusion constant was analytically 
derived for this model, showing that the divergence of the 
relaxation time at a finite density found in simulations is 
only apparent |l4| . The growth of t re \ (p) with p is in fact 
extremely abrupt, loglogt rc i(p) cx (1 — p)^ 1 , with an as- 
sociated coherence length loglog£(p) = c/(l — p), where 
c is a numerical constant. The extremely fast increase 
of t Te i can easily be numerically confused with a true di- 
vergence at a finite density. This more difficult problem 
represents thus a highly selective test for the methods 
described in this paper. We now show that finite size 
scaling allows one to distinguish between a true and an 
apparent divergence of the relaxation time in this model. 

We use again a continuous time Monte Carlo algorithm 
to study a model of linear size N = 24 for densities 
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FIG. 3: The parameter B(T, L) as a function of the re- 
duced variable \og(£(p)/L) with I — (p c — p)~ a (top) or 
I — expexp(c/(l — p)) (bottom). A much better collapse 
is obtained in the bottom part, with one less free parameter. 
Systematic deviations are visible in the top part. Both are 
for the kinetically constraint lattice glass of Ref. |4|. 



three arguments thus discriminate both possibilities, and 
we numerically confirm the absence of a dynamical arrest 
at finite density in this model. 

To summarize, we have shown that finite size scaling 
techniques can be used to study the glass transition. This 
follows from the identification of local two-time observ- 
ables F(r,t) as relevant 'order parameters' with interest- 
ing fluctuations and correlations. We have defined two 
quantities, B(T, L) and x(T, L), whose scaling behaviour 
allows one to locate the glass transition in a much more 
accurate manner than an extrapolation of simpler quan- 
tities such as the relaxation time. The method proposed 
here straightforwardly applies to any glass-forming sys- 
tems, including off-lattice models. More generally, our re- 
sults show that the glass transition has more in common 
with continuous phase transitions than usually believed 
and raise the hope that more concepts and methods bor- 
rowed from continuous phase transitions studies can also 
be used. 

I thank JP Garrahan and W Kob who provided me 
with their continuous time algorithms 0,0 , and G Biroli 
for communicating results before publication |l4j . Over 
the last 14 months, daily conversations on this subject 
with JP Garrahan were very useful. This work is sup- 
ported by a European Marie Curie Fellowship No HPMF- 
CT-2002-01927, CNRS (France) and Worcester College 
Oxford. Numerical results were obtained on Oswell at 
the Oxford Supercomputing Center, Oxford University. 



p E [0.75,0.87], covering almost 8 decades of relaxation 
times. We use F(i, t) = p(i, t)p(i, 0) as a dynamical quan- 
tity, where p(i,t) is the density at site i and time t. The 
averages are performed in subsystems of linear sizes L, 
using 10 4 independent initial conditions at each density. 
The correlator Q is well described by a simple exponen- 
tial form, C(r,t Te i) — e~ r / i ( p \ The latter measurements 
are also a check against spurious finite size effects, since 
they show that N 3> £(p) for the densities studied here. 

We show our results for the determination of the crit- 
ical density in Fig. |3 where the parameters B(p,L) are 
represented as a function of two possible scaling variables 
l/L. The top figure assumes a power law divergence of 
the coherence length t(p) w (p c — p)~ a at a finite density 
p c < 1, while the bottom figure makes use of the results 
of Ref. ^ij , with a divergence at p c — 1 only. The best 
collapses obtained in the two cases are shown, using the 
values p c = 0.882, a = 0.59, and c = 0.167. The second 
scaling is clearly superior, the first one having system- 
atic deviations either at large or small sizes. The figure 
shown here is a compromise between those. Moreover, 
the second scaling has one less free parameter. The same 
conclusions are drawn from the susceptibility ©. These 
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